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Abstract 

We consider the Blume-Capel model with the quantum tunneling between 
the excited states. We find a magnetically ordered phase transition induced 
by quantum fluctuation in a model. The model has no phase transition in 
the corresponding classical case. Usually, quantum fluctuation breaks ordered 
phase as in the case of the transverse field Ising model. However, in present 
case, an ordered phase is induced by quantum fiuctuation. Moreover, we 
find a phase transition between a quantum paramagnetic phase and a classi- 
cal diamagnetic phase at zero temperature. We study the properties of the 
phase transition by using a mean field approximation (MFA), and then, by a 
quantum Monte Carlo method to confirm the result of the MFA. 
Recently, in the study of functional material, various properties originating from the nontriv- 
ial energy structure of excited states have attracted much interest. In particular, when some 
cooperative interaction exists among the excited states, the system exhibits various phases 
which are incompatible to the ground state, and successive phase transitions occur among 
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them. For example, the effect of interaction on excitation has been studied with regard to 
mixed- valance biferrocenium salt, [1] and also the photo- induced magnetic properties. [2] 

The effect of the excited state is expressed by the Blume-Capel model with the variable 
Sf — 0, ±1 as the simplest model. [3] Here the state S^ — denotes the ground state and 
the states of Sf — ±1 denote the excited states with the magnetization or the polarization. 
In the Blume-Capel model, we assume a cooperative interaction among the excited states 
which could induces ordered states even in the cases where the state S- — has a lower local 
energy because of the kind of crystal field, etc. Even if the ordered state is not realized, the 
interaction still causes a metastable state of the magnetized phase (or the polarized phase), 
which could cause interesting photo-induced switching phenomena. [2] 

In the previous work, we introduced quantum tunneling between the states of — ±1 
on the Blume-Capel model in order to explain the phase transition by isotope replacement 
induced in SrTiOs. [4] Although the TiOe cluster of SrTiOa has the six excited states, we 
considered only the case of two excited states, that is, the Blume-Capel model, because 
the structure of the phase diagram does not change with the number of excited states. We 
also considered the tunneling only between the excited states because we assumed that the 
tunneling rate between the excited states was higher than the tunneling rate between ground 
state and the excited states in SrTiOs. 

In this paper, we study the two different effects of the quantum tunneling on this model. 
The first one is the effect as the quantum fluctuation. That is, as we know in the transverse 
Ising model, [5] the quantum fluctuation suppresses the ordered phase, and sometimes it 
destroys the order to cause the quantum disordered phase. We see this effect in the present 
model, too. The second one is the effect of the energy gain due to the zero temperature 
fluctuation. That is, the mixing of S'f = ±1 states reduces the ground state energy although 
it also reduces the total amount of magnetization. In the case of the ground state being in 
the state Sf — 0, this energy gain may cause the magnetized state to have a lower energy. 
In this case, an ordered phase appears in the region where an ordered state would not exist 
when the quantum effect did not exist. 
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In this paper we emphasize this second case as a new type of quantum effect on the 
magnetic ordered phase and investigate the phase transition. First, we consider the mean 
field approximation (MFA) and next, we perform a quantum Monte Carlo simulation to 
confirm the result of the MFA. 

The Hamiltonian of this model is given by 



(1) 



Here, is the z-component of the S — 1 spin operator and Sf is the operator which denotes 
the quantum tunneling between the excited states: 
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First, we treat this model by using the MFA. We replace by (S^) + Si and neglect the 
higher terms of 6. Then, the Hamiltonian is given by 



n = Y.i-Jmzst + A{sif + r^f ) - - W 
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where 



"^MFA 



r 





V 



(4) 



-Jzm + A - H 


r Jzm + A + H 

and z denotes the numbers of the nearest neighbor sites and H is the symmetry breaking 
field, m — (S^) stands for the order parameter. In this approximation scheme, the on-site 
quantum fluctuation is taken into exactly. Diagonalizing this Hamiltonian (4) , 
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the free energy is obtained as 
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FM¥A.{m) = -/cfiTlnTrexp (-/JT^mfa) - -Jzrri^ 



= -A;BTln(l + 2exp(-/3A) cosli(/3e)) - -Jzm% (6) 
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where e = y ( J^rm + i?)^ + F^. 

The self-consistent equation for m is 

_ 2{Jzm + //) exp(-^A) sinh(/?e) 
~ e(l + 2 exp(-/3A) cosh(/3e)) ' 

We solve this self-consistent equation numerically and we obtain the phase diagram as shown 
in Fig. 1. In the case of large A (A > j2;/2), a ferromagnetic phase appears. This ferromag- 
netic phase does not exist in the classical model, and is induced by quantum fluctuation. 
We investigate properties of the ground state in order to clarify the nature of the quantum 
fluctuation induced ordered phase. From the self- consistent equation for F < Jz, possible 
solutions for m are 



m = 0,±Wl-(i^y, (8) 



JzJ 

at T ^ and H ^ 0. For F > Jz, the (free) energy has only one minima at m = 0, 
where the quantum fluctuation destroys the magnetic order. For T < Jz, the energies of the 
system for the solutions (8) are 

E^O (m = 0) 

(9) 



We show the ground state phase diagram in Fig. 2. When A is small in the region F < Jz, 

the magnetic ordered state of m 7^ gives the lower energy states than that of m = 0. When 
we increase A, the ground state has a quantum phase transition at 

A=i^^fl+f^VV (10) 
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which is depicted by a solid curve. We notice that the critical value of A increases with the 
quantum fluctuation. In the classical case, although the ordered phase does not appear for 



^Jz < A. This phase boundary is indicated by horizontal dotted hne. Therefore, we find 
that the ordered phase is induced by the quantum fluctuation in the region ^Jz < A < Jz 
which is hatched. 

The quantum ferromagnetic state |Q.F.) is a mixed state of the up | t) and down | I) 
states: 

|Q.F.)=p|T) + ?|i), (11) 

where -\- (f — 1. For F < Jz, 

and = = I for F > Jz. That is, by the exchange interaction J, the mixing ratio of 
the up state and the down state is not the same. As the result, the quantum ferromagnetic 
phase appears. This type of mixed state appears also in the transverse field Ising model. 
However in the present case, the interesting thing is that the quantum fluctuation induces 
a magnetic order phase from the diamagnetic state where 5^ = is dominant. 

We also notice that at the quantum paramagnetic phase is different from the diamagnetic 
phase, although the magnetization is zero in both phases. Therefore, there is a quantum 
phase transition between these two phases at zero temperature. Although magnetization is 
always zero and does not show a jump, the susceptibility shows a jump at this quantum 
phase transition as shown in Fig. 3. 

In order to realize the present phase diagram, it is important that the tunneling between 
the excited states is larger than the tunnehng between the local ground state — Q and the 
local excited states m = ±1. For the comparison, we consider a model where the tunneling 
between the ground state and the excited states is given by the Hamiltonian 

n^-jY. sts] + A Y^istf + Tj2s:, (13) 

where Sf is the x component of 5" = 1 spin. As depicted in Fig. 4, the ordered phase is 
monotonically reduced with F and ordered phase induced by the quantum fluctuation does 
not appear. 
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Next, we study properties of the model (1) in the square lattice by a quantum Monte 
Carlo method with the Suzuki-Trotter decomposition. In the present model, non-diagonal 
interaction exists only at each site. Thus, the model is simply transformed to the classical 
cubic model. 

The Suzuki-Trotter decomposition transforms the Hamiltonian into a three-dimensional 
classical model TCeS- 



(14) 



TieS is given with a variable a = —1, 0, and 1 as 
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Here, N is the Trotter number and we adopt the periodic boundary condition for the Trotter 
direction ajv+i = We performed Monte Carlo simulation by using the cluster heat bath 
algorithm. [4,7] In Fig. 5, we present the phase diagram in the coordinate (F, T). In Fig. 
6, we present the ground state phase diagram in the coordinate (F, A). These results are 
quahtatively same with the results obtained by using the MFA in the previous section. Thus, 
we confirm qualitative features of the result obtained by the MFA. 

We found the ordered phase which is induced by the quantum phase transition. We 
think that this model is one of the prototype models which show the quantum fiuctuation 
induced phase transition. 

The numerical calculations were partially performed using the facilities of the Super- 
computer Center, ISSP, University of Tokyo. The present work is partially supported by a 
Grant-in- Aid from the Ministry of Education, Culture, Sports, Science and Technology. 
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FIG. 1. Phase diagram of the Blume-Capel model with the quantum tunnehng between the 
excited states in the coordinate (F, T). We choose the number of the nearest neighbor sites as 
z = 4. 
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FIG. 2. Ground state phase diagram in the coordinate (F, A). The solid line depict the first 
order transition and the dashed line depict the second order transition. We choose the number of 
the nearest neighbor sites as z = 4. 
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(a) (b) 

FIG. 3. (a) r dependence of the susceptibility at ksT/J = 0.0 and A = 5.0. (b) Temperature 
dependence of the susceptibihty at A = 5.0 and various T. The transition field is F = 5.0. 
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FIG. 4. Phase diagram of the Blume-Capel model with the quantum tunneling between the 
ground states and the excited states in the coordinate (F, T). We choose the number of the nearest 
neighbor sites as z = 4. 
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FIG. 5. Phase diagram in the coordinate (F, T) by using the quantum Monte Carlo simulation. 




FIG. 6. Ground state phase diagram in the coordinate (F, A) by using the quantum Monte 
Carlo simulation. The solid line depict the first order transition and the dashed line depict the 
second order transition. 
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